clear global;clc;
% 定义的各种系数
TriangleCoes;

% 使用高斯积分估计,
tef= @(para) testFTrans(para,matToXY,vec1);
% 1, 3, 4, 71, 72
ans_gauss=delta_e*TriangleGauss(71,tef)

% 使用 matlab integral2 积分
tef2=@(x,y) testFTrans2(x,y, matToXY, vec1);
ymax = @(x) 1 - x;
% Jacobi 因子是 2 delta_e
ans_matlab = 2*delta_e*integral2(tef2,0,1,0,ymax)


%+++++++++++++++++++++++++++++++++++++++++++
% 测试函数 f(x,y); 例如多项式
function F=testF(para)
x=para(1);
y=para(2);
F=x^3+y^3+6*x^2*y+6*x*y^2;
end

%  经过坐标变换
function F=testFTrans(para,mat,vec)
paraTmp=(mat*para'+vec)';
F=testF(paraTmp);
end

%  经过坐标变换
function F=testFTrans2(x,y,mat,vec)
[lenx,leny]=size(x);
F=zeros(lenx,leny);
for ni=1:lenx
    for nj=1:leny
        para=[x(ni,nj) y(ni,nj)];
        paraTmp=(mat*para'+vec)';
        F(ni,nj)=testF(paraTmp);
    end
end
end